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ABSTRACT 

A  computer  simulation  was  conducted  to  investi¬ 
gate  the  effects  of  both  linear  and  parabolic  tooth 
profile  modification  on  the  dynamic  response  of  low- 
contact-ratio  spur  gears.  The  effect  of  the  total 
amount  of  modification  and  the  length  of  the  modifica¬ 
tion  zone  were  studied  at  various  loads  and  speeds  to 
find  the  optimal  profile  modification  for  minimal 
dynamic  loading. 

Design  charts  consisting  of  normalized  maximum 
dynamic  load  curves  were  generated  for  gear  systems 
operated  at  various  loads  and  with  different  tooth 
profile  modification.  An  optimum  profile  modifica¬ 
tion  can  be  determined  from  these  design  charts  to 
minimize  the  dynamic  loads  of  spur  gear  systems. 


Ln  normalized  length  of  tooth  profile  modification 
zone  defined  such  that  Ln  =  1.0  is  the  length 
from  tooth  tip  to  HPSTC,  measured  along  the  line 
of  contact. 

R(j  base  radius,  mm  ( in . ) 

T(_  torque  on  load,  N-m  (in. -lb) 

T^  torque  on  motor,  N-m  (in. -lb) 

Tfi  torque  on  gear  1,  N-m  (in. -lb) 

Tf2  torque  on  gear  2,  N-m  (in. -lb) 

Wn  normalized  total  transmitted  load 


NOMENCLATURE 

Cg  damping  coefficient  of  gear  tooth  mesh,  N-sec 
(lb-sec) 

Cs  damping  coefficient  of  shaft,  N-m-sec 
( i n  .-lb-sec) 

J|_  polar  moment  of  inertia  of  load,  kg-m2 
(in. -lb-sec2) 

polar  moment  of  inertia  of  motor,  kg-m2 
( in. -lb-sec2) 

Jl  polar  moment  of  inertia  of  gear  1,  kg-m2 
( i n .-lb-sec2 ) 

polar  moment  of  inertia  of  gear  2,  kg-m2 
( i n . 1 b-sec2) 

kj  dynamic  factor 

kg  stiffness  of  gear  tooth,  N/m  (lb/in.) 
ks  stiffness  of  shaft,  N-m/rad  (in. -lb/rad) 


0  angular  di splacement,  rad 

9  angular  velocity,  rad/sec 

8  angular  acceleration,  rad/sec2 

A  amount  of  profile  modification  (thickness  of 

material  removed  from  tip  of  involute  gear 
tooth),  defined  such  that  A  =  1.0  is  the  mini¬ 
mum  amount  of  tip  relief  recommended  by  Welbourn, 

Vim 

INTRODUCTION 

One  of  the  major  concerns  in  the  design  of  power 
transmission  gears  is  the  reduction  of  gear  dynamic 
load.  Research  on  gear  noise  and  vibration  has 
revealed  that  the  basic  mechanism  of  noise  generated 
from  gearing  is  gear  box  vibration  excited  by  the 
dynamic  load.  Vibration  is  transmitted  through  shafts 
and  bearings  to  noise-radiating  surfaces  on  the  exte¬ 
rior  of  the  gear  box.  Dynamic  load  creates  cyclic 
bending  stresses  in  tooth  roots  which  can  lead  to 
fatigue  failure  as  well  as  cyclic  subsurface  stresses 
which  can  cause  tooth  surface  failure  by  pitting  and 
scoring.  The  life  and  reliability  of  a  gear  transmis¬ 
sion  is  reduced  by  high  dynamic  load.  Minimizing  gear 
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dynamic  load  will  decrease  gear  noise,  increase  effi¬ 
ciency,  improve  pitting  fatigue  life,  and  help  prevent 
gear  tooth  fracture  (1-5). 

Modifying  gear  tooth  profile  is  a  widely  used 
practice  to  reduce  dynamic  load  for  improved  perform¬ 
ance  of  a  spur  gear  transmission.  Current  practice  in 
gear  design  is  to  modify  the  tooth  profile  based  on 
the  maximum  applied  torque,  also  called  design  torque. 
When  a  modified  gear  system  is  operated  at  off-design 
torque,  its  dynamic  load  may  become  significant. 

Research  efforts  have  been  conducted  in  this  area 
for  many  years,  yet  there  is  a  lack  of  systematic  work 
leading  to  in-depth  understanding  of  how  tooth  profile 
modifications  affect  the  dynamic  response  of  a  spur 
gear  transmission  (1,4, 6-9). 

If  the  center  of  the  driven  gear  is  held  fixed 
and  a  torque  is  applied  at  the  center  of  the  driving 
gear,  the  teeth  in  contact  and  the  bodies  of  both 
gears  will  deform.  This  condition  yields  an  angular 
displacement  of  the  center  of  the  driving  gear  rela¬ 
tive  to  the  fixed  frame  of  reference  at  the  center  of 
the  driven  gear.  The  relative  angular  displacement  of 
the  gears  can  be  converted  to  a  linear  displacement 
along  the  line  of  action.  The  total  relative  displace¬ 
ment  of  the  driving  gear  with  respect  to  the  driven 
gear  along  the  line  of  action  is  defined  as  the  static 
transmission  error. 

This  paper  discusses  a  computer  simulation  study 
in  which  the  total  amount  of  tootn  profile  modifica¬ 
tion  and  the  length  of  the  modification  zone  were  sys¬ 
tematically  varied  to  determine  their  effect  on  the 
static  transmission  error  and  dynamic  loading  of  spur 
gears.  Both  linear  and  parabolic  modifications  were 
studied.  Their  individual  influence  and  relative  sig¬ 
nificance  on  gear  dynamic  load  are  compared  and 
di scussed. 

A  gear  set  which  operates  at  a  constant  design 
torque  can  be  optimally  modified  to  minimize  its 
dynamic  response.  For  gear  systems  that  are  to  be 
operated  at  off-design  torque  or  under  variable  load¬ 
ing  conditions,  design  charts  describing  the  gear 
dynamic  load  response  due  to  different  profile  modifi¬ 
cations  are  presented.  The  optimum  length  and  amount 
of  tooth  profile  modification  may  be  determined  from 
these  design  charts. 

THEORETICAL  ANALYSIS 

The  theoretical  model  assumes  that  a  simple  spur 
gear  transmission,  which  consists  of  a  driving  and  a 
driven  gears,  two  connecting  shafts,  a  motor,  and  a 
load,  can  be  treated  as  a  lumped-mass  vibration  system 
(Fig.  1)  (10,11).  The  motion  of  the  system  is  expressed 
by  the  following  set  of  differential  equations, 

Jm9M  +  C s i < ©M  -  9] )  +  Ksi(9m  -  9;)  -  TM 
J]9]  +  Cs](9]  -  9m )  +  K5i(9]  -  9^)  +  Cg(t) 

CRbl©l  -  Rb292 1  ♦  Kg<t)  CRbl (Rbl9l  -  Rb292>3  »  Tf 1 <  t) 

J2S2  +  Cs2<92  ‘  ®l)  +  Ks2(92  -  ®L>  +  cg<*> 

CRb292  -  Rp 1 ©1 ]  +  Kg(t>  CRb2(Rb2e2  -  Rb19l)]  ■  Tf2(t> 

JL9L  *  (9[_  -  *2>  ♦  Ks2(9L  -  92)  -  -tl 

Where  9m,  9j ,  92 ,  and  9^  represent  the  rota¬ 
tions  of  the  motor,  the  gears,  and  the  load;  Jm,  Ji. 

J2,  and  represent  the  mass  moments  of  Inertia 
of  the  motor,  the  gears,  and  the  load;  Csi.  Co.  and 
Cg ( t )  are  damping  coefficient  of  the  shafts  and  the 
gears;  Ksi,  k,2,  and  Kg(t)  are  stiffnesses  of  the 


shafts  and  the  meshing  teeth.  Tm,  T ,  Tf](t),  and 
Tf2(t)  are  motor  and  load  torques  and  frictional 
torques  on  the  gears;  Rpi  and  Rq2  are  base  circle 
radii  of  the  gears;  t  is  time;  and  the  dots  over  sym¬ 
bols  indicate  time  differentiation. 

In  developing  the  above  equations  several  simpli¬ 
fying  assumptions  were  employed.  The  dynamic  process 
is  defined  in  the  rotating  plane  of  the  gear  pair,  and 
the  contact  between  gear  teeth  is  assumed  to  be  along 
the  theoretical  line  of  action.  Damping  due  to  lubri¬ 
cation,  etc.  is  expressed  as  a  constant  damping  factor 
(ratio  of  the  damping  coefficient  to  the  critical  damp¬ 
ing  coefficient.)  From  gear  research  literature,  typ¬ 
ical  damping  factors  of  0.10  and  0.005  respectively 
were  chosen  for  the  tooth  mesh  and  and  for  the  con¬ 
necting  shafts  (12  to  14). 

The  stiffnesses  and  mass  moments  of  inertia  of 
the  system  components  were  found  from  the  fundamental 
mechanics  of  materials  principles.  The  equations  of 
motion  contain  the  excitation  term  due  to  periodic 
variation  of  the  mesh  stiffness  and  due  to  errors 
(such  as  spacing  or  profile  errors).  The  meshing 
stiffness  is  a  function  of  the  mesh  point  along  the 
1 .ne  of  action.  Detailed  analyses  of  the  tooth  mesh¬ 
ing  stiffness,  shared  tooth  load,  and  static  transmis¬ 
sion  error  of  the  meshing  gear  pair  were  presented  in 
previous  studies  (9,10). 

Figure  2  presents  a  flowchart  of  the  generalized 
computational  procedure  for  the  solution  of  the  govern¬ 
ing  differential  equations.  The  equations  were  linear¬ 
ized  by  dividing  the  mesh  period  into  small  intervals. 

A  constant  input  torque  tm  was  assumed.  The  output 
torque  TL  was  considered  to  fluctuate  as  a  result  of 
time-varying  stiffness,  friction,  and  damping  in  the 
mesh. 

To  start  the  solution  iteration  process,  initial 
values  of  the  angular  displacements  were  obtained  by 
preloading  the  input  shaft  with  the  nominal  torque 
carried  by  the  system.  Initial  values  of  the  angular 
speed  were  taken  from  the  nominal  system  operati  ig 
speed. 

The  iterative  procedure  was  as  follows,  the  cal¬ 
culated  values  of  the  angular  displacement  and  speed 
after  one  mesh  period  were  compared  with  the  assumed 
initial  values.  Unless  the  differences  between  them 
were  smaller  than  a  preset  tolerance,  the  procedure 
was  repeated  using  the  average  of  the  initial  and  cal¬ 
culated  values  as  new  initial  conditions.  More  com¬ 
plete  descriptions  of  this  method  may  be  found  in 
Refs.  9  and  10  and  similar  work  appears  in  Refs.  4,  6 
and  7 . 

The  analysis  was  applied  to  a  sample  set  of  gears 
as  specified  in  Tabie  I.  These  are  identical  low- 
contact-ratio  spur  gears  with  solid  gear  bodies.  The 
number  of  teeth  is  28  and  the  module  is  3.18  mm.  Face 
width  is  25.4  mm  with  a  design  load  of  350  000  N/m 
(2000  Ib/in).  The  gear  mesh  theoretical  contact  ratio 
is  1.64.  A  typical  gear  tooth  showing  both  the  unmod¬ 
ified  (true  involute)  and  modified  profiles  is  illus¬ 
trated  in  Fig.  3(a).  A  sample  profile  modification 
chart  is  shown  in  Fig.  3(b).  On  the  chart,  a  straight 
line  represents  a  linear  tooth  profile  modification  and 
a  parabolic  line  represents  a  parabolic  modification. 

In  this  study,  the  same  amount  and  the  same 
length  of  profile  modifications  were  applied  to  the 
tooth  tip  of  both  pinion  and  gear.  The  minimum  amount 
of  conventional  tip  relief  was  chosen  as  a  reference 
value  to  normalize  the  amount  of  profile  modification. 
Hence,  for  the  minimum  amount  of  conventional  tip 
relief,  A  =  1.00.  Welbourn  stated  that  the  minimum 
tip  relief  should  be  equal  to  twice  the  maximum  spac¬ 
ing  error  plus  the  combined  tooth  deflection  evaluated 
at  the  highest  point  of  single  tooth  contact  (HPSTC) 
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(15).  The  length  of  profile  modification  is  designated 
Ln.  The  distance  along  the  tooth  profile  from  tooth 
tip  to  the  HPSTC  is  defined  to  be  of  unit  length.  The 
values  of  A  and  L„  can  be  varied  arbitrarily  to 
obtain  any  desired  combinations 

Figure  3(b)  shows  examples  of  linear  and  parabolic 
profile  modifications.  In  both  cases  the  amount  of 
modification  A,  .  1.00,  and  the  modification  length 
Ln  «  1.00.  Although  the  length  of  modification  is 
shown  as  a  vertical  distance  parallel  to  the  tooth 
axis  in  Fig.  3(a),  it  is  actually  defined  in  terms  of 
the  gear  roll  angle  as  specified  in  Fig.  3(b). 

RESULTS  AND  DISCUSSION 

Figures  4  and  5  show  the  comparison  of  the  static 
transmission  errors  s^d  shaded  tooth  ’oad:  for  oripc.1'  ■ 
fied  ge;rs  and  those  with  linear  and  parabolic  tooth 
profile  modifications.  The  normalized  modification 
length  Ln  was  set  at  1.0,  which  means  the  tip  relief 
extended  from  tooth  tip  to  the  HPSTC  location.  The 
modification  amount  varied  from  A  =  0.50  to  A  =  1.25 
at  an  increment  of  0.25.  When  the  amount  of  profile 
modification  was  less  than  or  equal  to  the  minimum 
conventional  tip  relief,  (A  £  *  1.00),  the  length  of 
single  and  double  contact  zones  as  shown  on  the  static 
transmission  error  graphs  were  not  changed  and  the  con¬ 
tact  ratio  remained  at  1.64.  When  an  excessive  modi¬ 
fication  amount  (for  example,  A  =  1.25)  was  applied  on 
the  tooth  profile  the  zone  of  double  tooth  contact 
shortened  and  gear  contact  ratio  was  reduced  (to 
approximately  1.53  for  this  case). 

The  principal  excitation  for  gear  system  vibra¬ 
tion  is  the  unsteady  component  of  the  re’ative  angular 
motion  of  meshing  gears  due  to  the  variation  of  static 
transmission  error.  (The  steady  part  of  transmission 
error  which  is  due  to  gear  body  "windup"  does  not 
cause  excitation.)  The  main  purpose  of  profile  modi¬ 
fication  is  to  minimize  this  variation.  A  comparison 
of  the  conventional  tip  relief  curves  (A  »  1.0, 
ln  =  1.0)  in  the  static  transmission  error  plots  of 
Figs.  4  and  5  shows  that  the  linear  profile  modifica¬ 
tion  curve  is  smoother  than  the  one  with  parabolic 
modification.  This  indicates  that  if  the  conventional 
amount  and  length  of  tip  relief  is  used,  a  spur  gear 
system  with  linear  profile  modification  is  expected  to 
provide  a  smoother  dynamic  response  than  gears  with 
parabolic  modification. 

Figure  6  shows  a  speed  sweep  plot  of  the  dynamic 
load  factor  for  gears  with  no  tip  relief  (unmodified), 
and  gears  with  linear  and  parabolic  tooth  profile  mod¬ 
ifications  (conventional  modification  amount  and 
length;  A  =  1.0,  Ln  *  1.0).  The  dynamic  load  factor 
is  defined  as  the  ratio  of  maximum  dynamic  tooth  load 
during  contact  to  static  tooth  lead.  The  primary  res¬ 
onance  for  these  cases  occurs  near  fundamental  system 
natural  frequency  11,280  rpm.  A  Jacobi  iteration 
technique  was  (16)  used  to  determine  the  system  natu¬ 
ral  frequencies.  The  peak  value  for  the  unmodified 
case  was  about  2.18.  Peak  values  for  the  linear  and 
parabolic  cases  were  approximately  1.30,  and  1.40 
respectively.  As  above,  the  linear  tip  relief  yields 
the  smoothest  response. 

To  understand  the  detailed  effect  of  tooth  pro¬ 
file  modification  on  the  dynamic  behaviour  of  a  spur 
gear  transmission,  the  amount  (A)  and  length  (Ln)  were 
varied  systematically.  First,  the  effect  of  linear 
tooth  profile  modification  on  the  dynamics  of  the  sam¬ 
ple  gears  was  investigated.  Figure  7  shows  the  speed 
sweep  of  dynamic  load  factor  for  the  sample  gear  sys¬ 
tem  with  linear  tooth  profile  modification  running  at 
design  load.  The  normalized  length  was  ln  =  1.00  and 
the  amount  was  varied  from  A  =  0.75  to  A  -  1.25. 


The  dynamic  response  of  a  unmodified  gear  pair  is  also 
shown  for  comparison.  As  expected,  the  peak  dynamic 
load  factor  at  resonance  speed  is  minimum  at  A  =  1.0 
and  ln  =  1.0.  The  maximum  dynamic  effect  at 
A  =■  0.75  was  less  than  that  at  A  =  1.25.  This  result 
was  anticipated  in  Fig.  4(a)  where  there  is  less  vari¬ 
ation  in  the  static  transmission  error  curve  at 
A  =  0.75  than  at  A  »  1.25.  This  last  result  suggests 
that  there  is  a  greater  detrimental  effect  of  excess 
profile  modification  than  of  under  modification. 

Excess  profile  modification  reduces  the  contact  ratio 
which  increases  the  dynamic  load. 

Figure  3  shews  the  effect  of  varying  load  on  the 
dynamic  response  of  the  sample  gear  set  with  conven¬ 
tional  linear  tip  relief  (A  =  1.0,  Ln  =  1.0).  In 
Fig.  8(a),  the  applied  load  was  normalized  using  the 
design  load  ( 350  000  N/m)  as  me  reference  value. 

(Ie;  Wn  =  1.00  when  the  appled  load  equals  the  design  ; 

load.)  At  design  load  (Wn  «  1.00),  the  value  of  the 

peak  dynamic  factor  is  1.30.  This  is  the  minimum 

dynamic  factor  found.  As  the  applied  load  varies  from 

the  design  load,  the  maximum  dynamic  load  factor 

increases  from  this  value.  From  Fig.  8,  load  factor 

curves  are  shown  at  normalized  applied  load  values 

(Wn)  of  0.6,  0.8,  1.0,  and  1.2.  The  corresponding 

peak  values  of  dynamic  load  factor  are  approximately 
2.47,  1.82,  1.30,  and  1.45.  These  curves  also  show 
that'underload  (Wn  <1 .0)  produces  a  greater  dynamic 
load  factor  than  overload  (Wn  >1.0).  Finally,  the 
dynamic  load  factor  curve  of  an  unmodified  (involute) 
gear  pair  under  design  load  is  shown  for  comparison. 

The  peak  dynamic  load  value  for  unmodified  involute 
gears  is  2.18. 

To  obtain  a  more  realistic  feeling  of  the  actual 
dynamic  loading  on  the  gear  tooth  and  to  prevent  mis¬ 
leading  interpretation,  the  speed  sweep  curves  of 
Fig.  8(a)  were  replotted  to  show  the  actual  tooth  load 
in  Fig.  3(b).  The  smallest  peak  value  of  the  dynamic 
tooth  load  occurs  for  the  design  load  case  (Wn  =  1.0). 

Both  underload  and  overload  cases  show  higher  values 
of  peak  load.  From  the  curves,  maximum  dynamic  load 
values  found  are  518,700  N/m,  510,000  N/m,  455,000  N/m, 
and  609,000  N/m  for  Wn  =  0.60,  0.80,  1.00,  and  1.20 
respectively.  The  detrimental  effect  of  operating 
gears  at  a  load  substantially  lower  than  the  design 
load  and  at  the  resonant  speed  was  clearly  demonstrated. 

For  at  Wn  =  0.60,  the  peak  dynamic  tooth  load  »as 
actually  greater  than  that  at  Wn  =  0.80  and  Wn  =  1.00. 

Once  again,  the  curve  for  unmodified  involute  gears 
running  at  Wn  =  1.0  is  shown  for  comparison.  The 
benefit  of  gear  tooth  profile  modification  can  be  seen 
by  comparing  the  dynamic  tooth  loads  of  modified  and 
unmodified  gears. 

Similar  studies  were  performed  on  the  sample 
gears  with  parabolic  tooth  profile  modif ications .  The 
results  are  presented  in  Figs.  9  and  10.  Unlike  the 
linear  modification  case,  the  minimum  dynamic  response 
of  gears  with  parabolic  profile  modifications  with 
ln  =  I. 00  occured  at  A  -  1.25  instead  of  at  A  =  1.0. 

This  can  be  explained  by  comparing  the  static  trans¬ 
mission  error  curves  of  these  two  cases  in  Fig.  5(a). 

At  Ln  =  1.00,  the  error  curve  for  A  ■  1.25  is 
smoother  than  that  for  A  =  1.0.  This  means  gears 
with  parabolic  tooth  profile  modifications  require  a 
greater  amount  of  modification  than  gears  with  linear 
profile  modifications. 

Figure  10  shows  the  dynamic  response  curves  of 
gear  pairs  modified  with  paracolic  tip  relief  at 
A  -  1.00  and  Ln  =  1.0  for  various  applied  loads.  The 
curve  at  Wn  =  0.8  had  the  lowest  peak  value.  Contrary 
to  the  linear  case,  gears  vd th  parabolic  modifications 
run  more  smoothly  at  underload  than  at  design  load. 
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From  the  above  observation,  one  may  conclude  that 
for  conventional  amount  and  length  of  profile  modifi¬ 
cation  (A  «  1.00  and  Ln  »  1.00),  linear  profile  modi¬ 
fication  should  be  used  for  gears  which  will  operate 
at  and  above  design  load,  and  parabolic  profile  modi¬ 
fication  should  be  applied  to  gears  operating  below 
design  load,  to  minimize  dynamic  effect. 

The  various  effects  of  applied  load,  profile  mod¬ 
ification  length,  and  profile  modification  amount  on 
the  normalized  maximum  dynamic  load  of  spur  gears  with 
either  linear  or  parabolic  tooth  profile  modifications 
were  further  investigated.  The  noramlized  maximum 
dynamic  load  is  defined  as  the  product  of  the  maximum 
dynamic  load  factor  (MDLF)  and  the  normalized  total 
transmitted  load  (Wn).  This  normalized  magnitude  of 
the  maximum  dynamic  load  in  the  gear  mesh  provides 
better  comparison  of  gear  dynamics  at  different 
applied  loads.  Multiplying  this  normalized  value  by 
the  design  load  gives  the  actual  gear  dynamic  load. 

First,  a  constant  modification  length  of 
Ln  =  1.00  was  assumed,  and  three  different  modifica¬ 
tion  amounts  of  i  =  0.75,  1.00,  and  1.25  were 
applied  to  the  sample  gears.  In  Fig.  11  are  plotted 
curves  of  the  normalized  maximum  dynamic  load  over 
the  load  range  of  0.70  to  1.20  times  the  design  load 
(Wn).  For  the  linear  modification  case,  shown  in 
Fig.  11(a),  the  normalized  maximum  dynamic  load  reaches 
a  minimum  value  at  0.76  Wn  on  the  A  =  0.75  curve 
and  at  1.00  Wn  on  the  A  =  1.00  curve.  The  minimum 
of  the  A  =  1.25  curve  apparently  occurs  at  a  load 
greater  than  1.2  Wn  and  is  therefore  off  the  scale 
of  Fig.  11(a).  The  normalized  maximum  dynamic  load 
appears  to  be  more  sensitive  to  load  change  for  over¬ 
load  than  for  underload. 

Figure  11(b)  presents  the  dynamic  load  data  for 
the  parabolic  modification  case.  On  the  curves  for 
A  -  0.75  and  A  »  1.00,  the  minimum  dynamic  effect 
occurs  at  a  load  less  than  Wn  .  0.70  and  thus  off  the 
scale  of  Fig.  11(b).  On  the  curve  for  A  *  1.25,  the 
minimum  occurs  at  approximately  0.72  Wn. 

Comparing  the  curves  in  Figs.  11(a)  and  (b)  shows 
that  the  gears  with  parabolic  tip  relief  are  much  less 
sensitive  to  changes  in  the  amount  of  tip  relief  than 
gears  with  linear  tip  relief.  Therefore,  it  is 
expected  that  the  dynamics  of  parabolic  tip  relieved 
gears  would  be  less  affected  by  manufacturing  toler¬ 
ances  and  machining  errors.  In  addition,  the  normal¬ 
ized  maximum  dynamic  load  for  gears  with  parabolic 
relief  appears  to  be  generally  lower  than  for  gears 
with  linear  relief  over  the  load  range  of  Hn  «  0.7  to 
1.2  (underload  to  overload).  This  means  that  parabolic 
tip  relief  is  clearly  a  better  choice  than  linear  tip 
relief  for  gears  that  must  operate  over  a  wide  range 
of  loads. 

The  effect  of  different  amounts  of  profile  modi¬ 
fication  on  the  normalized  maximum  dynamic  load  of 
gears,  at  various  applied  loads  in  the  range  of 
Wn  »  0.7  to  Wn  =  1.2,  is  shown  in  Fig.  12.  As  in  the 
previous  figure,  the  length  of  the  modification  zone 
was  held  constant  at  ln  »  1.00.  Figure  12(a)  shows 
the  curves  for  gears  with  linear  modifications,  and 
Fig.  12(b)  for  those  with  parabolic  modifications. 

The  optimum  amount  of  profile  modification  for  gears 
operating  at  either  a  single  load  or  over  a  range  of 
loads  can  be  estimated  from  the  minimum  points  on 
these  curves.  For  the  linear  modification  case, 

A  =  i.00  is  optimum  for  gears  operating  at  the  design 
load  (constant  Wn  =  1.0).  If  the  gears  operate  over 
a  range  of  loads,  the  optimum  amount  of  modification 
is  found  from  the  Intersection  of  the  curves  for  the 
highest  and  lowest  loads  of  the  range.  Therefore,  for 
loads  ranging  from  Wn  .  0.7  to  Wn  =  1.0,  the  optimum 
modification  occurs  at  A  -  0.92  which  corresponds  to 


the  intersection  of  the  Wn  =  0.7  and  Wn  =  1.0  curves 
in  Fig.  12(a).  Likewise,  A  =  1.18  is  optimum  for 
gears  that  operate  from  Wn  =  0.7  to  Wn  =  1.2.  For 
the  parabolic  modification  case,  it  appears  that 
A  ■  1.25  is  the  optimum  amount  for  gears  operating 
from  Wn  =  0.7  to  either  Wn  =  1.0  or  Wn  =  1.2.  As 
noted  above  in  the  discussion  for  Fig.  11,  the  dynamic 
response  of  parabol ical ly  modified  gears  is  less 
affected  by  the  changes  in  the  amount  of  profile  modi¬ 
fication  than  are  gears  with  linear  modification. 

Finally,  the  effect  of  length  of  tooth  profile 
modification  on  spur  gear  dynamic  response  was  inves¬ 
tigated  and  is  shown  in  Fig.  13.  The  modification 
amount  was  held  constant  at  A  =  1.00.  The  length  of 
modification  zone  was  varied  from  Ln  .  0.50  to  1.30 
and  maximum  dynamic  load  curves  were  generated  for 
several  values  of  applied  load  (Wn).  The  minimum 
dynamic  response  for  gears  with  linear  tooth  profile 
modification  occured  at  Ln  =  0.67,  0.78,  1.00  respec¬ 
tively  for  Wn  <=  0.70,  0.80,  and  1.00,  see  Fig.  13(a). 
Since  gears  seldom  operate  at  a  constant  load  in  their 
daily  operation  a  method  must  be  found  to  choose  pro¬ 
file  modification  specifications  for  the  selected 
design  load  range.  For  the  load  range  of  0.70  to  1.00 
of  design  load  (0.7  <  Wn  <1.0),  an  optimum  length  for 
linear  tooth  profile  modification  is  Ln  *  0.90.  This 
value  is  obtained  from  the  intersection  point  of  the 
Wn  =  0.70  and  Wn  *  1.00  curves  from  the  normalized 
maximum  dynamic  load  curves  in  Fig.  13(a).  Any  modi¬ 
fication  length  other  than  this  would  yield  less 
desirable  higher  dynamic  effect  under  this  range  of 
loads. 

A  similar  study  for  parabolic  tooth  profile  modi¬ 
fication  is  shown  in  Fig.  13(b).  The  applied  load  was 
varied  from  0.70  to  1.20  of  design  load.  (This  is  a 
wider  load  range  than  used  for  the  linear  case  above, 
since  we  have  shown  that  gears  with  parabolic  modifi¬ 
cations  are  suitable  for  a  wider  load  range.)  An 
optimum  length  of  modification  for  minimum  dynamic 
response  for  gears  operating  over  a  range  of  loads  may 
be  determined  from  this  figure.  For  example:  At  con¬ 
stant  design  load,  (Wn  •  1.0),  the  optimum  length  of 
modification  is  Ln  =  1.30.  For  overload  (Wn  >  1.0), 
the  curves  suggest  that  the  optimum  length  will  be 
greater  than  1.30  (thus  extending  beyond  the  pitch 
point).  In  this  study,  modifications  extending  beyond 
the  pitch  point  were  not  considered.  As  another 
example,  if  the  operating  load  range  is  Wn  =  0.70  to 
Wn  -  1.00  (underload  to  design  load),  the  optimum 
length  is  approximately  Ln  =  1.28  (found  at  the 
intersection  of  the  Wn  =  0.70  and  Wn  =  1.00  curves). 
Finally,  for  a  wider  load  range  of  Wn  =  0.70  to 
Wf,  =  1.20  (underload  to  overload),  the  length  of  modi¬ 
fication  is  chosen  to  be  1.30  (since  this  study  does 
not  consider  modification  extending  beyond  the  pitch 
point).  In  general,  a  longer  (than  1.0)  length  of 
modification  zone  seems  to  be  preferred  for  parabolic 
tooth  profile  modification  since  it  yields  lower 
dynamic  load. 

A  comparison  of  figures  12  and  13  shows  that  the 
modification  length  (Ln)  has  a  greater  Impact  on  the 
maximum  dynamic  load  factor  than  does  the  amount  of 
modification  (A).  Therefore  the  length  of  modifica¬ 
tion  should  be  controlled  as  closely  as  possible. 
Nevertheless,  due  to  machining  errors  and  allowable 
tolerance  it  is  not  practical  to  manufacture  tooth 
profile  modifications  exactly  as  specified  by  the 
theory.  In  reality,  a  modified  tooth  profile  devi¬ 
ates  somewhat  from  the  ideal  specification.  As  dis¬ 
cussed  earlier,  parabolic  profile  modification 
appears  to  be  less  sensitive  to  manufacturing  vari¬ 
ance  and  is  therefore  preferred  to  linear  profile 
modification. 
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As  an  example  of  designing  the  optimum  parabolic 
tooth  profile  for  a  spur  gear  transmission  operating 
at  a  range  of  loads,  consider  a  gearset  wnich  operates 
over  the  load  range  between  Wn  =  0.7  and  Wn  =  1.2. 
Since  the  dynamic  load  is  more  sensitive  to  the  length 
of  modification  < Ln )  than  to  the  amount  (A),  ln  is 
chosen  first.  From  figure  13(b)  the  optimum  length  is 
1.30.  With  the  length  Ln  fixed  at  this  value,  the 
optimum  amount  of  profile  modification  can  be  found  by 
varying  A  over  a  suitable  range  as  shown  in  figure  14. 
This  figure  shows  dynamic  load  curves  at  applied  loads 
<Wn)  of  0.7,  l.C,  and  1.2  for  gears  with  modification 
length  Ln  =  1.30,  and  modification  amount  varying 
from  A  =  0.75  to  4  =  1.50.  The  optimum  amount  of 
profile  modification  is  found  to  be  A  =  1.18.  This 
is  the  intersection  point  of  the  Wn  =  0.7  and  the 
Wn  »  1.2  curves.  For  this  example,  the  worst  case 
(highest  value)  of  normalized  maximum  dynamic  load 
will  be  1.40.  This  is  the  load  corresponding  to  the 
extremes  of  the  range  of  applied  load  (at  Wn  =  0.70 
and  at  Wn  =  1.20). 

CONCLUSIONS 

A  computer  simulation  was  conducted  to  investi¬ 
gate  the  effects  of  both  linear  and  parabolic  tooth 
profile  modifications  on  the  dynamic  response  of  low- 
contart-ratio  spur  gears.  The  effects  of  the  total 
amount  of  modification  and  the  length  of  the  modifica¬ 
tion  zone  were  studied  at  various  loads  and  speeds  to 
find  optimal  (low  dynamic  response)  specifications  for 
profile  modification.  The  following  conclusions  were 
obtained : 

1.  The  amount  and  type  of  tooth  profile  modifica¬ 
tions  have  a  significant  effect  on  the  dynamic  per¬ 
formance  of  spur  gear  systems. 

2.  Parabolic  tooth  profile  modification  is  gener¬ 
ally  preferred  for  low  dynamic  response  in  gears  which 
operate  over  a  range  of  loading  conditions.  These 
gears  are  less  sensitive  to  changes  in  applied  load, 
amount  of  modification  and  length  of  modification  than 
are  gears  with  linear  profile  modifications. 

3.  Gears  with  parabolic  profile  modifications 
require  a  slightly  longer  length  of  modification  zone 
than  gears  with  linear  profile  modifications.  The 
modification  zone  may  extend  beyond  the  highest  point 
of  single  tooth  contact. 

4.  Gears  which  operate  at  a  nearly  constant  load 
at  design  load  to  moderate  overload  will  perform  more 
quietly  (with  less  dynamic  effect)  if  linear  profile 
modification  is  used. 

5.  For  gears  with  linear  profile  modification, 
excess  modification  has  a  greater  detrimental  effect 
on  dynamic  loads  than  under  modification,  and  under¬ 
load  causes  higher  dynamic  effect  than  overload. 

6.  Over  a  range  considered  in  this  report,  the 
length  of  modification  has  a  greater  effect  on  the 
dynamic  response  for  both  linear  and  parabolic  pro¬ 
file  modifications  than  does  the  total  amount  of 
modi f i cation . 
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TABLE  I.  -  GEAR  DATA 


Gear  tooth  .  .  .  Standard  involute  fu 
Module,  mm  (diametrial  pitch,  in.-1) 

11-depth  tooth 

.  .  .  3.18  (8) 

Pressure  angle,  deg  . 

. 20 

Number  of  teeth  . 

. 28 

Face  width,  mm  (in. ) . 

.  .  25.4  (1 .0) 

Design  load,  N/m  (lb/in.)  . 

350  000  (2000) 

Theoretical  contact  ratio  . 

. 1.64 

( 


0 


L 


(a)  A  SIMPLE  GEAR  TRANSMISSION. 


(b)  SYMBOLIC  MODEL. 

FIGURE  1.  -  COMPUTER  MODEL  OF  SPUR  GFAR  SYSTEM. 
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GEOMETRY  OF  SYSTEM  COMPONENTS  AND 
CONDITIONS  OF  SYSTEM  OPERATION 


CALCULATION  OF  TOOTH  PROFILE, 
MESH  STIFFNESS.  INERTIAS. 
DAMPING.  AND  FRICTION 


CALCULATION  OF 
INITIAL  CONDITIONS 


CALCULATION  OF  DYNAMIC  CONDITION 
FOR  ONE  PERIOD  PER  MESH 


ASSUMPTION  OF  NEW 
INITIAL  CONDITION 


CONVERGENCE  OF 
DYNAMIC  CONDITION? 


CALCULATION  OF  DYNAMIC 
LOADS  AND  STRESSES 


OUTPUT  OF  RESULTS 


FIGURE  2.  -  FLOW  CHART  OF  COMPUTATIONAL  PROCEDURE. 
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NORMALIZED  AMOUNT  OF  TOOTH 
PROFILE  MODIFICATION,  A 


V 

\ 

LENGTH  OF  PROFILE  \ 


MODIFICATION,  Ln  \ 


TRUE  INVOLUTE  TOOTH  PROFILE 


^AMOUNT  OF  PROFILE 
MODIFICATION,  A 


■MODIFIED  PROFILE 


HIGHEST  POINT  OF 


SINGLE-TOOTH 
,  l  CONTACT,  HPSTC 

\  l 

'  '-PITCH  POINT 


LOWEST  POINT  OF  SINGLE-\ 
TOOTH  CONTACT,  LPSTC 


(a)  GEAR  TOOTH  WITH  MODIFIED  TOOTH  PROFILE. 


IS  18  20  22  24  2G  28  30 

ROLL  ANGLE,  DEG 

fb)  SAMPLE  TOOTH  PROFILE  MODIFICATION  CHART. 


FIGURE  3.  -  EXAMPLE  OF  MODIFIED  GEAR  TOOTH. 
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TOOTH  LOAD,  LB/ IN.  TRANSMISSION  ERROR.  IN. 
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(b)  SHARED  TOOTH  LOAD. 


FIGURE  4.  -  STATIC  TRANSMISSION  ERROR  AND  SHARED 
TOOTH  LOAD  FOR  GEAR  PAIRS  WITH  LINEAR  TOOTH  PRO 
FILE  MODIFICATIONS.  FULL  DESIGN  LOAD;  LENGTH  OF 
MODIFICATION,  Ln  =  1.00. 
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TOOTH  LOAD,  lb/in.  TRANSMISSION 
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(a)  STATIC  TRANSMISSION  ERROR. 
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(b)  SHARED  TOOTH  LOAD. 

FIGURE  5.  -  STATIC  TRANSMISSION  ERROR  AND  SHARED 
TOOTH  LOAD  FOR  GEAR  PAIRS  WITH  PARABOLIC  TOOTH 
PROFILE  MODIFICATIONS.  FULL  DESIGN  LOAD; 
LENGTH  OF  MODIFICATION,  Ln  =  1.00. 
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DYNAMIC  LOAD  FACTOR  DYNAMIC  LOAD  FACTOR 


FIGURE  6.  -  DYNAMIC  LOAD  FACTOR  OF  SPUR  GEAR  PAIRS 
UNDER  DESIGN  LOAD  WITH  TRUE  INVOLUTE  TOOTH  PRO¬ 
FILE,  LINEAR  PROFILE  MODIFICATION,  AND  PARABOLIC 
PROFILE  MODIFICATION. 


FIGURE  7.  -  EFFECT  OF  VARYING  AMOUNT  OF  LINEAR  TOOTH 
PROFILE  MODIFICATION  ON  DYNAMIC  LOAD  FACTOR  OF 
SPUR  GEAR  PAIR.  FULL  DESIGN  LOAD;  LENGTH  OF  MODI¬ 
FICATION,  Eft  =  1.00 
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(a)  DYNAMIC  LOAD  FACTOR. 


0  3000  6000  9000  12  000  15  000 

SPEED,  RPM 

(b)  TOTAL  DYNAMIC  TOOTH  LOAD. 


FIGURE  8.  -  EFFECT  OF  VARYING  APPLIED  LOAD  ON  DYNAMIC 
LOAD  FACTOR  AND  TOTAL  DYNAMIC  LOAD  OF  SPUR  GEAR 
PAIR.  CONVENTIONAL  LINEAR  TIP  RELIEF;  LENGTH  OF 
PROFILE  MODIFICATION,  Ln  =  1.0;  AMOUNT  OF  PROFILE 
MODIFICATION.  A  =  1.0. 
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FIGURE  9.  -  EFFECT  OF  VARYING  AMOUNT  OF  PARABOLIC 
TOOTH  PROFILE  MODIFICATION  ON  DYNAMIC  LOAD  FAC¬ 
TOR  OF  SPUR  GEAR  PAIR.  FULL  DESIGN  LOAD;  LENGTH 
OF  MODIFICATION,  Ln  =  1.00. 


SPEED,  RPM 

FIGURE  10.  -  EFFECT  OF  VARYING  APPLIED  LOAD  ON  DYNAMIC 
LOAD  FACTOR  OF  A  SPUR  GEAR  PAIR.  PARABOLIC  TIP 
RELIEF;  LENGTH  OF  PROFILE  MODIFICATION,  Ln  =  1.0; 
AMOUNT  OF  PROFILE  MODIFICATION,  A  =  1.0.  (UNMODI¬ 
FIED  INVOLUTE  CASE  IS  ALSO  SHOWN  FOR  COMPARISON.) 
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(b)  PARABOLIC  PROFILE  MODIFICATION. 


FIGURE  11.  -  EFFECT  OF  APPLIED  LOAD  ON  NORMALIZED  MAXIMUM 
DYNAMIC  LOAD  OF  SAMPLE  GEARS  AT  VARIOUS  MODIFICATION 
AMOUNT.  LENGTH  OF  PROFILE  MODIFICATION,  Ln  =  1.00. 
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(b)  PARABOLIC  PROFILE  MODIFICATION. 


FIGURE  12.  -  EFFECT  OF  PROFILE  MODIFICATION  AMOUNT 
ON  NORMALIZED  MAXIMUM  DYNAMIC  LOAD  OF  SAMPLE 
GEARS  AT  VARIOUS  APPLIED  LOADS.  LENGTH  OF  PRO¬ 
FILE  MODIFICATION,  Ln  =  1.00. 
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FIGURE  13.  -  EFFECT  OF  LENGTH  OF  PROFILE  MODIFICATION 
ON  NORMALIZED  MAXIMUM  DYNAMIC  LOAD  OF  SAMPLE  GEARS 
AT  VARIOUS  APPLIED  LOADS.  AMOUNT  OF  PROFILE  MODI¬ 
FICATION,  A  =  1.00 
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FIGURE  14.  -  OPTIMUM  PARABOLIC  PROFILE  MODIFICATION 
FOR  SAMPLE  GEARS  OVER  RANGE  OF  APPLIED  LOADS. 
LENGTH  OF  PROFILE  MODIFICATION.  Ln  =  1.30. 
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